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\ Abstract 

The Vlasov-Poisson-Boltzmann System governs the time evolution of the distribu- 
tion function for the dilute charged particles in the presence of a self-consistent electric 
potential force through the Poisson equation. In this paper, we are concerned with the 
J> ■ rate of convergence of solutions to equilibrium for this system over M 3 . It is shown that 

£N| [ the electric field which is indeed responsible for the lowest-order part in the energy 

space reduces the speed of convergence and hence the dispersion of this system over 
the full space is slower than that of the Boltzmann equation without forces, where 
the exact difference between both power indices in the algebraic rates of convergence 
^1 \ is 1/4. For the proof, in the linearized case with a given non-homogeneous source, 

■ Fourier analysis is employed to obtain time-decay properties of the solution operator. 

f— ^ In the nonlinear case, the combination of the linearized results and the nonlinear en- 

ergy estimates with the help of the proper Lyapunov-type inequalities leads to the 
optimal time-decay rate of perturbed solutions under some conditions on initial data. 
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1 Introduction 

The Vlasov-Poisson-Boltzmann (called VPB in the sequel for simplicity) system is a phys- 
ical model describing the time evolution of dilute charged particles (e.g., electrons) under 
a given external magnetic field |2H [1] . The VPB system for one-species of particles in the 
whole space M 3 reads 

dtf + Z- V,./ + V x .$ • V 6 f = Q(f, /), (1.1) 
A x $ = [ fd£-p(x), (1.2) 

with initial data 

/(0,x,0=/o(*,0- (°) 
Here, the unknown / = f(t,x,£) is a non-negative function standing for the number den- 
sity of gas particles which have position x = (x±,X2, £3) £ K 3 and velocity £ = (£1, £2, £3) £ 
M 3 at time t > 0. Q is the bilinear collision operator for the hard-sphere model defined by 



Qif, 9)= / if 9* -/<?*)!(£ -Z*)-u\dwd&, 

Jr 3 xS 2 

f = fit,x,0, f' = fit,x,C'), g* = g(t,x,&), g'*=git,x,Q, 
£' = £-[(£-£*) & = + 

The potential function $ = x) generating the self-consistent electric field in (jl.lj) is 
coupled with /(t,x,£) through the Poisson equation (|1.2|) . p(a;) denotes the stationary 
background density satisfying 

/?(x) — ► /Oqo as |x| — > 00, 

for a positive constant state p^ > 0. 

The existence of stationary solutions to the system (|1.1|) - (|1.2|) and the nonlinear stabil- 
ity of solutions to the Cauchy problem (|l.l|) - ([1.3p near the stationary state were obtained 
in [8], and the corresponding results have been recalled in the appendix. In this paper, 
we are concerned with the rate of convergence of solutions towards the stationary states. 
Since the background density does not produce any essential difficulty, for simplicity it is 
supposed throughout this paper that 

p(x) = Poo = i, xem 3 . (1.4) 

In this case, the VPB system (jl.ip - (|1.2|) has a stationary solution (/*, <!>*) with 

/*=M, $* = 0, 



where 
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is a normalized global Maxwellian in three dimensions which has zero bulk velocity and 
unit density and temperature. One of the main results of this paper is stated as follows. 
Notations and norms will be explained at the end of this section. 



Theorem 1.1. Let N > A and iu(£) = (1 + |£| 2 ) 1/2 . Assume that f > and 

/o-M 



M 



H N nLl,nz 1 



(1.5) 



is sufficiently small. Let f > be the solution obtained in Proposition \5.2\ to the Cauchy 
problem (jl.ip . (jl.2p and (| 1 .3|) under the assumption (|1.4p . Then, f enjoys the estimate 
with algebraic rate of convergence: 



fit) - M 


< c 


/o-M 






y/M 



H N r\Z! 



(!+*)■ 



(1.6) 



Furthermore, under the following additional conditions on /o, / also enjoys some estimates 
with extra rates of convergence: 



Case 1. If 



[ (/ (x,O-M)de = 
holds for any x G ]R 3 , t/ien one /ias 



/(*) - M 


< c 


/o-M 









(!+*)• 



(1.7) 



H N nz 1 



Case 2. Fix any < e < 3/4 and suppose that 

/o-M 



M 



(1.9) 



is sufficiently small. Then one has 



f(t)-T,U(^J(t)) ej M 



M 



ii 1 



+|| (eo, f{t) - M) + l|V..(e„ /(i)}||^ 

3=1 



< c 



/o-M 



M 



(! + *)■ 



(1.10) 



where {ej}1- =0 is the orthonormal set in L 2 (]R 3 ; M<i£) defined by 



eo = 1, e.j = Cj (1 < j < 3), e 4 



V6 ' 



4 



R.-J. Duan and R.M. Strain 



Remark 1.1. The rates of convergence in (jl.6p and (|1.8p are optimal under the corre- 
sponding assumptions in the sense that they coincide with those rates given in Theorem \3.1\ 
at the level of linearization. Similarly, the rate of convergence in (jl.lOp is almost optimal. 
To our knowledge, these results are the first ones in the study of rate of convergence for 
the nonlinear VPB system in IR 3 . Of course, it is an interesting problem to improve the 
almost optimal rate in (jl.lOp to the optimal rate. 

Remark 1.2. In \32j on the same issue, some slower algebraic rates of convergence of 
solutions in the space L^°(L|) were obtained on the basis of the pure energy estimates and 
time-decay of some ordinary differential equation. Here, under the assumptions of (|TT 
from the Sobolev inequality, (jl.lOp shows that 



f(t) - M 



< C 



/o-M 



M 



(i + *)-!- 



where L^(H^~ 2 ) = L°°(IR 3 ; H N ~ 2 (M.^; wd£)). The above inequality implies the almost 
optimal rate of convergence of solutions. The main reason why one can here obtain the 
optimal or almost optimal rates is that we make full use of the time-decay properties of 
solutions to the linearized system with nonhomogeneous sources. This will be carried out 
in Section 0, where we shall also state another main result Theorem \3.1\ in this paper. 



Remark 1.3. As shown in 129], for the Boltzmann equation without external forces, the 
power index in (jl.6|) takes the value 3/4. Hence, the dispersion of the VPB system inM 3 is 
slower than that of the Boltzmann equation. This is essentially caused by the self-induced 
potential force. Furthermore, if one decomposes f as the summation of three parts 



4 [ 4 

/ = (e , />e M + f) ej M +lf- f) ej M 

j'=i [ i=o 

then by comparing (jl.6p and (jl.lOp . one can find out that the effect of the self-induced 
potential force on rates of convergence only happens in the above second part which cor- 
responds to the projections of f along the momentum components ej (1 < j < 3) and the 
temperature component e±. This phenomenon is consistent with that recently obtained by 
\1 7f for the Navier-Stokes-Poisson system in M 3 . 

The time rate of convergence to equilibrium is an important topic in the mathematical 
theory of the physical world. As pointed out in [30], there exist general structures in 
which the interaction between a conservative part and a degenerate dissipative part lead 
to convergence to equilibrium, where this property was called hypocoercivity. Theorem ll.il 
indeed provides a concrete example of the hypocoercivity property for the nonlinear VPB 
system in the framework of perturbations. The key of the method to study hypocoercivity 
provided by this paper is to carefully capture the time-decay rates for the perturbed 
macroscopic system of equations with the hyperbolic-parabolic structure, which is in the 
same spirit of the Kawashima's work [16]. 
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There has been extensive investigations on the rate of convergence for the nonlinear 
Boltzmann equation or related spatially non-homogeneous kinetic equations with relax- 
ations. In what follows let us mention some of them. In the context of perturbed solutions, 
the first result was given by Ukai [28], where the spectral analysis was used to obtain the 
exponential rates for the Boltzmann equation with hard potentials on torus. The results in 
|28| were improved by Ukai- Yang [29] in order to consider existence of time-periodic states 
in the presence of time-periodic sources, which was later extended by Duan-Ukai- Yang- 
Zhao [10] to the case with time-periodic external forcing by using the energy-spectrum 
method. We also mention Glassey-Strauss [12] for the study of the essential spectra of 
the solution operator of the VPB system. Recently, Strain-Guo [27] developed a weighted 
energy method to get the exponential rate of convergence for the Boltzmann equation 
and Landau equation with soft potentials on the torus. Earlier but along the same line 
of research, Strain-Guo [26] developed a general theory of polynomial decay rates to any 
order in a unified framework and applied it to four kinetic equations, the Vlasov-Maxwell- 
Boltzmann System, the relativistic Landau-Maxwell System, the Boltzmann equation with 
cutoff soft-potentials and the Landau equation all on the torus. 

Another powerful tool is entropy method which has general applications in the exis- 
tence theory for nonlinear equations. By using this method as well as the elaborate analysis 
of functional inequalities, time-derivative estimates and interpolation, Desvillettes- Villain 
[2] obtained first the almost exponential rate of convergence of solutions to the Boltzmann 
equation on torus with soft potentials for large initial data under the additional regularity 
conditions that all the moments of / are uniformly bounded in time and / is bounded 
in all Sobolev spaces uniformly in time. See Villani [30] for extension and simplification 
of results in [2] still conditionally to smoothness bounds by further designing a new aux- 
iliary functional. Notice that [26] provided a simple proof of [2] for the unconditional 
perturbative regime. Recently, by finding some proper Lyapunov functional defined over 
the Hilbert space, Mouhot-Neumann [23] obtained the exponential rates of convergence 
for some kinetic models with general structures in the case of torus; see also [30] for the 
similar study. An extension of [23] to models with additional confining potential forces 
was given by Dolbeault-Mouhot-Schmeiser [5]. 

Besides those methods mentioned above for the study of rates of convergence, the 
method of Green's functions was also founded by Liu-Yu [20] to expose the pointwise 
large-time behavior of solutions to the Boltzmann equation in the full space ffi 3 . 

Here, we mention that if there is no collision in (|l.l[) - (|1.2p . that is to consider the 
Vlasov-Poisson system, then the so-called Landau damping comes out. This was recently 
studied by Mouhot- Villani [24] on torus, where it was shown that even though in the 
absence of kinetic relaxation, in the analytic regime the solution still converges weakly in 
some sense to certain large-time states determined by the initial data and the nonlinear 
system itself, for any interaction potential less singular than Coulomb. In the Coulomb 
case, they established Landau damping over exponentially long times. 

Finally, we also mention some of results on the existence theory of the VPB system and 
related kinetic equations: global existence of renormalized solutions with large initial data 
BUSES], global existence of classical solutions near Maxwellians [13 [02 EJJ , OH USES] 
and [H |8] , and global existence of solutions near vacuum [HI [9] . 
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The rest of this paper is organized as follows. In Section [2] we make some preparations 
to reformulate the Cauchy problem of the VPB system, make the macro-micro decom- 
position for both the solution and equations, obtain a system of equations describing the 
evolution of some macroscopic velocity moment functions for the later analysis in both the 
linear and nonlinear cases, and reduce Theorem 11.11 in an equivalent form to Proposition 
12.11 In Section [3] and Section 01 we obtain the time-decay rates of perturbed solutions 
under some conditions on initial data in the linear and nonlinear cases, respectively. Here, 
Theorem 13.11 which is another main result of this paper is applied together with some 
energy estimates to prove Proposition 12.11 Some free energy functionals and Lyapunov- 
type inequalities play a key role in the proof of Theorem 13. II and Proposition 12.11 Finally, 
we conclude this paper with an appendix in Section [5] by listing some results obtained in 
[8] about the existence of the stationary solution and its nonlinear stability for the VPB 
system (fTT]l -(fl~2l). 

Notations. Throughout this paper, C denotes some positive (generally large) constant 
and A denotes some positive (generally small) constant, where both C and A may take 
different values in different places. In addition, A ~ B means \\A < B < A2^4 for two 
generic constants Ai > and A2 > 0. For an integrable function g : M n — > K, its Fourier 
transform g = Tg is defined by 

r. n 

g(k) = Tg(k) = / e~ 2M g(x)dx, x ■ k =: V Xjkj, 

for k £ R n , where i = £ C is the imaginary unit. For two complex vectors a,b £ C", 

(a I b) = a ■ b denotes the dot product over the complex field, where b is the complex 
conjugate of b. For any integer m > 0, we use H m , H™, H™ to denote the usual Hilbert 
spaces H m (R% x R"), H m (W^), H m (R™), respectively, where L 2 , L 2 ,, L 2 are used for the 
case when m = 0. For a Banach space X, \\ ■ \\x denotes the corresponding norm, while || • || 
always denotes the norm || ■ \\ L 2 or || ■ \\ L 2 for simplicity. We also use || • \\[j m for the norm of 
the weighted Hilbert space H m (R™ x k|; w(Z)dxd£) and || • || L 2 for L 2 (M.fx w(£)dxd£). 
We use (•, ■) to denote the inner product over the Hilbert space L?, i.e. 

(g,h) = I gmmt, g,h€Ll 

JRri 

For q > 1, we also define 



Z q = q{Ll)=L*(B%-,L*(B»)), \\g\\z q =([ (f 



\g(x,0\ q dx) dd 



For the multiple indices a = (a\, • • • , a n ) and (3 = ((3\, • • • , j3 n ), as usual we denote 

°x °£ —° Xl °X n °^ % • 

The length of a is \a\ = a± + ■ ■ ■ + a n . For simplicity, we also use dj to denote d x . for each 
j = !,■■■ ,n. Generally, except in Section [3l we consider the case of n = 3 dimensions. 
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2 Macro-micro decomposition 

In this section, we shall make some preparations for the later analysis in both the linear 
and nonlinear cases, and moreover reduce Theorem ll.il to Proposition 12. ll in the equivalent 
form. Firstly, one can reformulate the Cauchy problem (jl.ip . (jl.2p and (jl.3p as follows. 
Set the perturbation u = u(t,x,^) by 



/ = M + VMu. (2.1) 
Then u and $ satisfy the perturbed system: 

d t u + £ • V x u + V a .$ • V 5 u - -£ • V x $u - V x .$ • £\/M = Lu + T(u, tt), (2.2) 

Aa;$ = / VMud£, (2.3) 

with given initial data 

u(0,x,O=u o (x>S) = ^=¥ i , (2-4) 

V JVC 

where Lit and r(u, u) are denoted by 

1 



Lu 



Q(M, a/Mu) + Qt^Mu, M)l , (2.5) 



T(u, u) = -^Q(VMu, VMu). (2.6) 



It is well-known that for the linearized collision operator L, one has 

(Ku)(£) = [ (-VMu* + y/M^u' + VM'u'*) |(£ - £*) • ti}\y/M^dud€* 

where z/(£) is called the collision frequency and K is a self-adjoint compact operator on 
L 2 (M|) with a real symmetric integral kernel -?T(£,£*). The null space of the operator L 
is the five dimensional space spanned by the collision invariants 

AA = J fCerL = span|v / M; &VM,i = 1,2,3; |£| 2 \/m} . 

From the Boltzmann's H-theorem, the linearized collision operator L is non-positive and 
moreover, — L is locally coercive in the sense that there is a constant A > such that 

-[ uLud£>\ [ z/(£) ({I - ~P}u) 2 d£, VueD(L), 

JR 3 Jr 3 
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where for fixed (t,x), P denotes the projection operator from L\ to N and D(L) is the 
domain of L given by 

D(L) = {uGL 2 | y^ueLiy 
In addition, for the hard sphere model, v satisfies 

"(0 ~ (i + l£l 2 )* =™(0- 

This property will be used throughout this paper. 

Given any u(t,x,£), we define the projection operator, Pu, as 



Pu = <j a u (t,x) + ^2b^(t,x)^ + c u (t,xm 2 \ VM. 



(2.7) 



Since P is a projector, it holds that 

/ (i,e,iei 2 )VM{i-p« = o, 

Jw? 

i.e. {I — P}u is orthogonal to jV, which together with the form (|2.7p of P imply 

«" = \l (5-|£lVMud£, (2. 
2 Jr3 

6 U = f tVMud£, 
c U = \ I (\e-3)VMud£, 

where b u = 62,63). Thus, u(i, can be uniquely decomposed into 
' n(t, x, = Pu {I - P}u, 
Pu = {a u (t, x) + b u (t, x)-C + c u (t, x)|£| 2 } VM G A/", 
{I-P}ueAA\ 



(2.9) 
(2.10) 



(2.11) 



where Pu is called the macroscopic component of u(t,x,£) with coefficients (a u ,b u ,c u ), 
and {I — P}u the microscopic component of u(t,x,£). For later use, one can further 
decompose Pu as 



Pu = P U © PlU, 

P u = (a u + 3c u )VM, a u + 3c u = [ 

JR3 



Mud£, 



(2.12) 



P 1 n = {6"-e + c ?t (|e| 2 -3)}VM 



where Po and Pi are the projectors corresponding to the hyperbolic and parabolic parts 
of the macroscopic (fluid) component, respectively. 
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In what follows, we shall apply the macro-micro decomposition (|2.1ip to the system 
(|2.2p - (|2.3p to deduce the macroscopic balance laws satisfied by (a u ,b u ,c u ). Firstly, by 
multiplying (11,10 by the collision invariants 1,£, |£| 2 , one can get the local conservation 
laws 



dt / M + V x • / Udi = 0, 

Jr 3 Jm 3 

cm + v x • / e ® cm - / m = o, 



ft/ Id'M + v*-/ I^M-v^-/ ?M = o. 

JR3 J K 3 J E 3 



Plugging / = M + V MPu + yM{I — P}u into the above equations as well as into the 
Poisson equation (|1.2|) gives 

d t (a u + 3c u ) + V x -b u = 0, (2.13) 
+ V^o" + 5c 11 ) + V x • (£ ® £>/M, {I - P}u) - V a $ = (a u + 3c tt )V*$, (2.14) 
8 t (3a u + 15c") + 5V X • 6" + V x • m^VM, {I - P}u) = b u • V x $, (2.15) 

and 

A x $ = a u + 3c u . (2.16) 
Here and hereafter we use the moment values of the normalized global Maxwellian M: 

(1,M) = 1, 

<|£/,M> = 1, <|£| 2 ,M)=3, 
(\C j \ 2 \U\ 2 ,M) = l, jjLm, 

<y 4 ,M> = 3, <|£|%| 2 ,M> = 5, <|£| 4 ,M> = 15, 
<|£|%| 2 ,M>=35, <|£| 6 ,M> = 105. 

Notice that (|2T3j) and (|2T5j) implies 

d t c u + \v x ■ b u + \v x ■ UtftiVM, {I - P}u) = 0, 
3 o 

which is more convenient to be used to replace the time derivative of c u in the proof. 
As in [7], we also have to consider the evolution of higher-order moments of {I — P}u: 

(^(Vm,{i-ph <|£|Vm,{i-p}u>, 

which have appeared in (|2.14j) and (|2.15p , respectively. The equation (|2.2p can be rewritten 
as 

d t Pu + i ■ V x Pu - V x $ • tVM = -d t {I - P}u + R + G, (2.17) 
where the linear term R and the nonlinear term G are denoted by 

R = _f • V X {I - P}u + L{I - P}u, (2.18) 
G = T(n, u) - • V{« + • V x &) u. (2.19) 
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One can use (|2.1ip to further write (|2.17p as 

d t a u VM + ^2{d t b] + d ja u - dj^VM + Y,i d tc u + 0j&"}|&| 2 VM 
j j 
+ + dmb^UVM + E d jC u \C\%VM 

j<m j 

= -dt{l - P}u + R + G. (2.20) 

Define the high-order moment functions A = (>lj m )3x3 and B = (B\, B2, -B3) by 

A jm {u) = ((^ m -l)VM,u), B j ( U ) = m 2 -5)£ j VM,u). (2.21) 

Applying Aj m (-) and Bj(-) to both sides of (12.2QH and using the conservation law of mass 
(|2.13p . one has 

dtiAjjdl - P}u) + 2c u ] + 28jb] = Ajj (R + G), (2.22) 
d t A jm ({l-'P}u)+d j b^ + d m b] = A jm (R + G), j^m, (2.23) 
d t Bj({I - P}u) + d jC u = Bj(R + G), (2.24) 

where (I2.23P also holds for j > m since it is symmetric for (j, m) due to the symmetry of 
Aj m . The main observation which initially came from [13] and later [7] is that for fixed 
m, from (|2.22|) and (|2.23p . one can deduce 



^djA jm {{l - P}u) + -d m A mm ({I - P}u) 



^xbm dmdmbr, 



= \Y. d rnA ]3 {R + G)-J2 d,A im (R + G). (2.25) 

jj^m j 

Remark 2.1. It should be pointed out that the proof of (|5.1ip in Lemma l5.il for the 

macroscopic dissipation is only based on similar moment equations in presence of external 
forcing corresponding to (l2TT3l) - (l2TT5]l . (|2T6jK (|232>(|2T24l ) and (|235j) . The derivation of 
these moment equations was inspired by [13$ and firstly given by ffjjj and in the case of 
the Boltzmann equation. They play a key role in the choice of the solution space without 
time derivatives. Actually, the method of excluding time derivatives in 0, [?| is much 
useful since it will also be applied in Section of this paper to the study of the linearized 
system in order to obtain the time-decay property of solutions. 

Next, as mentioned at the beginning of this section, we claim that Theorem 11.11 is 
indeed implied by the following proposition in terms of perturbation u. Hence, the rest of 
this paper is devoted to the proof of this proposition. 

Proposition 2.1. Let N > 4 and = (1 + |£| 2 ) 1/2 . Assume that f = M + VMu > 
and H^oll/f^nL 2 nz 1 ^ s sufficiently small. Let f = M+\/Mu > be the solution obtained in 
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Proposition \5.2\ to the Cauchy problem (|2.2p . (| 2 . 3 1) and (|2.4|) under the assumption (jl.4|) . 
Then, u enjoys the estimate with algebraic decay rate in time: 

\\u{t)\\ HN < C \\u Q \\ HNnZl (1 + t)~l (2.26) 

Furthermore, under the following additional conditions on uq, u also enjoys some estimates 
with extra decay rates in time: 

Case 1. If 

P u = (2.27) 

holds for any i£l 3 , then one has 

\\u(t)\\ HN <C\\u \\ HNnZi (l + t)-l (2.28) 

Case 2. Suppose that \\uq\\ h n C}Zi is further sufficiently small. Then, for any < e < 3/4, 
there is rj > depending only on e such that whenever 

Ikollijjvnzi ^ r ?' 

one has 

||{I - Px}u(t)\\ HS + || V*Piu(t) < C \\uq\\hs^l (! + *r !+e - ( 2 - 29 ) 

Proof of Theorem ll.lt By the definition (12. ip for perturbation u, (11. 6p and (jl.8h are 
equivalent with (|2.26p and (|2.28j) . respectively, and conditions (|1.5p and (jl.9p coincide 
with those in Proposition I2.lt an d the condition (jl.7p is equivalent with (|2.27p since by 

P n = / VMu d£VM= [ (f -M)d^VM 

Jr 3 Jw 3 

holds for any (x,£) £ M 3 x M 3 . Finally, (jl.lOp is implied by (|2.29j) . This can be seen from 

||{I - Pi}u(t)||^ ~ ||{I - P}u(t)\\ H N + ||Pow(t)||i| W) , 
and further from 

4 

P n = (e , / - M)e Vm, P lU = ^(ej, f)e 3 VM, 

3=1 

and 

4 

{I - P}u = u- ^2(ej, VMu)ejVM 

3=0 

4 

J2(ej,f -M) ej VM 



M 

3=0 

^4 



/-Ej=o(ei,/)eiM 



/M 

L4 



Here, the orthonormal set {ej}| =0 is defined by (|1.11|) . Therefore, Theorem II .11 is proved. 
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Remark 2.2. Let us explain that conditions of Proposition POl make sure those of Propo- 
sition [57B are satisfied. Actually, this follows from the fact that 

8{uq) < C\\u \\ 2 H N nZl , 

where £{uq) is defined in (j5.3l) . It suffices to verify 

HVa-Aj^ouof < C\\u \\ 2 L2nZi . (2.30) 

In fact, it follows from the definition of Po and an interpolation inequality that 

HV^A^Pouoll 2 ~ ||V^o|| 2 < CUpollillpollli, 

X X 

where <J>o an d Po are defined by 



Furthermore, one has 



Then, (|2~3U|) follows. 



$0 = -T~T7 * /°o> Po= VMu d£. 
4ir\x\ J R 3 



\pohi < Ikoll, WpoWli < Ikobj- 



3 The linearized system 

In this section, let us consider the Cauchy problem of the linearized system with a non- 
homogeneous source: 



d t u + £ • V x u - • £VM = Lu + h, 

A x ^= I VMud£, t > 0,x £ IR n ,£ £ K n , 
^ «|t=o = ^o, a; £ K n , 

where h = h(t,x,t;) and uq = Uq{x,£) are given, the spatial dimension n > 1 is supposed 
to be arbitrary in order to see how it enters into the time-decay rate at the level of 
linearization, and for simplicity we still use M to denote the normalized n-dimensional 
Maxwellian 

M = - e"l«l 2 / 2 

(2vr)-/ 2 

Formally, the solution to the Cauchy problem (|3.1[) can be written as the mild form 



u (t) = e tB u + f e (t - s ^ B h(s)ds, (3.2) 
J 

where e iB denotes the solution operator to the Cauchy problem of the linearized equation 
without source corresponding to (13. ip for h = 0. The goal of this section is to show that 
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e has the proposed algebraic decay properties as time tends to infinity. The idea of 
proofs is to make energy estimates for pointwise time t and frequency variable k which 
corresponds to the spatial variable x. To the end, for 1 < q < 2 and integer m, set the 
index (T q>m of the time-decay rate by 



n 



>q,m 



1 1 



77? 



As mentioned at the end of Remark 11.21 another main result of this paper is stated in the 
following 



Theorem 3.1. Let 1 < q < 2 and n > 1. Let P and Pq be defined in (|3.6p . 

(i) For any a, a' with a' < a, and for any uq satisfying d^uo G L? and d^Uo G Z q , one 
has 



and 



\d«e tB u \\ + H^V^A-^oe^noll < C(l + t)~^-' (||^'u ||z 9 + ||<9>o||), (3-3) 

||c£e' B {I - P }no|| + Wd^V^V^il - P }u || 

< C (l + t)-^(\\d«'{I - P }u \\z q + \\d%{I - P }«o||), (3-4) 

for t > with m = \a — a'\, where C is a positive constant depending only on n,m,q. 

(ii) Similarly, for any a, a' with a' < a, and for any h satisfying ^(£) -1 / 2 <9"/i(/j) G L 2 and 
vfa-W&jf h{t) G Z q for t > 0, one has 

rt 2 rt 2 

d% I e^- s)B {I-P}h(s)ds + d*V x A~ l P e {t - s)B {I-P}h(s)ds 
Jo Jo 

(3.5) 



<C (l + t- s )- 



(\\^ 2 d«'{I - P}h(s)f Zq + ^-^{1 - P}h(s)f)ds, 
for t > with m= \a — a'\, where C is a positive constant depending only on n,m,q. 

Remark 3.1. In the case of the Boltzmann equation jl28f or the Navier-Stokes-Poisson 
system \17^ , some similar algebraic time-decay estimates in Theorem \3. 1\ were obtained on 
the basis of the spectral analysis of the solution semigroup. However, so far there has not 
been known results applying the direct spectral analysis as in [28] to the study of the VPB 
system. Here, the proof of Theorem \3.1\ that we shall show can provide a robust method 
to get the time-decay estimates in I? space for not only the VPB system but also some 
other kinetic models such as the Boltzmann equation and Landau equation, relativistic or 
non-relativistic, which is a future research goal in the general framework as in 130^ . 



Remark 3.2. In the case when the spatial domain is a torus T n , Glassey- Strauss 111)/ and 
Mouhot- Neumann f23\j obtained the exponential time-decay rate 



|e tB {I- P}u \\ x < Ce~ xt \ 



mo be, 
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where X = L 2 (T n x R n ) in JTT^ and X = ^(TP 



in 



Hf. See also JM GE ^ for 



the other kinetic models on torus. Actually, the proof of Theorem \3.1\ can be modified in 
a simple way so that the above inequality with exponential rates is recovered. We shall 
further explain this point at the end of this section; see Theorem \3.1\ 



To prove Theorem 13. 1( let u(t), formally defined by (j3.2j) . be the solution to the 
linearized non-homogeneous Cauchy problem (|3.ip . where uo(x, £) and h(t, x, £) with Ph 
I ) are given. We decompose u as 



( u(t,x,0 = Pu® {I-P}u, 
Pu = P u © P lU = {a u + b u -C + c u \Z\ 2 }VM, 
P u = (a u + nc u )VM, 
{P lU =[b u -C + c u (\C\ 2 -n)]VM, 
where a u ,b u , c u are the macro moment functions of u given by 

a u = \! [(n + 2)-|e| 2 ]^M^, 

b u = 



(3.6) 



£VMud£, 

? = 7T I (\H\ 2 -n)VMudZ, 
2n J R 3 



which are the n dimensional generalizations of (12. 8ft . (|2.9j) . (12.10j) when n = 3. Notice that 
although a u , b u and c u also depend on the spatial dimension, we used the same notations 
as in Section [2] for simplicity. Then, from the same procedure as in Section [21 one has the 
macroscopic balance laws satisfied by a u ,b u ,c u : 



d t (a u + nc u )+V x -b u = 0, 

d t b] + dj(a u + nc u ) + 2d jC u + d m A jm ({I - P}u) - fy® = 0, 



d t c u + -V x -b u + ^-Y] djBjdl - P} 



n 

a u + nc v 



0. 



(3.7) 
(3.8) 

(3.9) 
(3.10) 



and 



d t [A j:j ({I - P}u) + 2c u ] + 2djb] = AjjiR + h), 
d t A jm ({I - P}u) + djb^ + d m b] = Aj m (R + h), j + 



dtBj({I - P}u) + d jC u = Bj(R + h), 
for 1 < j, m < n. Here, the velocity moment functions Aj m (-) and Bj(-) are given by 



(3.H) 
(3.12) 
(3.13) 



A jm (u) = ((ZjU - l)VM,u), Bj{u) = ~{n + 2)&VM, 
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where these moment functions correspond to (|2.2ip for n = 3, and we again used the same 
notations as in Section [2] for simplicity. R has the same form with (|2.18p . given by 



(3.14) 



R = -£ . V X {I - P}u + L{I - P}u. 



Notice that the source term h does not appear in the first n + 2 equations (|3.7| )- (l3.8j) due 
to P/i = 0. Furthermore, similar to derive (|3.15[) from (|2.22[) and (|2.23|) . it follows from 
(IBTTTD and ([312]) that 



-dt 



J2djA jm ({I - P}u) + -d m A mm ({I - P}v 

3 

= lj2 9 m A lj(R + h)~J2 d 3 A 3m{R + h) 



A h u - 8 8 h u 



(3.15) 



Lemma 3.1. There is a free energy functional Ef ree (u(t,k)) which is local in the time 
and frequency and takes the form of 



E free (u(t, k)) 

E (e ifw Ajmi{1 - p ^ + 2ifw Ammi{I ~ P}S) 



Kl 



lk 3 A. (fX Dl^\ , 1 lkm A (fX Dl ^ I b U 



+«iX)(5 i ({I-P}S)| T 



ik 



■3 c u 



+E(^irrw (aS+ncG) 

for some constant k\ > 0, such that one has 

\k\ 2 

d t ReE free (u(t, k)) + 



(3.16) 



lu 1 2 i I u 1 2 



u i _, n 1 2 



A|a u + nc 



^c^iki-ph^ + ii^/^i-p}^ 



(3.17) 



for any t > and A; S M n . 



Proof. We shall make estimates on b u , c u and a u + nc" individually and then take the 
proper linear combination to deduce the desired free energy inequality (|3. 17j) . Firstly, 
notice that 

Ta u = a Tu , Th u = b^ u , Fc u = c Tu , 
and likewise for the high-order moment functions Aj m (-) and Bj(-). 
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Estimate on b u . We claim that for any < 8\ < 1, it holds that 



9 t Re Yl 2 ik 3 A jm({i- ~ P}S) + hk m A mm {{l - P}i 



< 5i(l + |fc| 2 )|a s + nc s | 2 + 5i|£f |c 



|2i mi2 



c 



[I-P}S|| 2 +|| I ,-V2 { I_p}fe||i 



In fact, the Fourier transform of (|3.15|) gives 



(3.18) 



-a 



ikjA jm ({I - P}u) + -ik m A mm {{\ - P}5 

= o X] ik mAjj(R + h) -y^JkjA jm (R + h). 



Taking further the complex inner product with 6^ gives 



+ (l*| 2 + Ol&*l 2 



d t ]J2 ikjA jm {{I - P}u) + Uk rn A rnm ({I - P} 
l - Y ik m Ajj{R + h)-^2 ikjA jm (R + h) | b q m 

i¥=rn j J 

+ ikjA jm ({I - P}u) + Uk m A mm ({I - P}2) | -d t bl j = Ji + J 2 . (3.19) 



C 



ii is bounded by 

h < 5i|Af lCl 2 + j Y.(\ A im(R)\ 2 + \A jm (h)?)- 

For I2, one can use the Fourier transforms of (|3.8|) and (|3.10D : 

d t bf + iAy K + (" + 2)c s ] + ^ ifc m A im ({I - P}u) - %8 = 0, 



4 $ = a u + rac u 



to estimate it as 



C 



(3.20) 
(3.21) 



h <h{l + \k\ 2 W + nc u \ 2 + 5^1 V| 2 + -(1 + \k\ 2 ) £ |A im ({I - P} 



On the other hand, notice from (13. 14ft that 

R = -£ • fc{I - P}u + L{I - P}u 
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which implies 

\A jm (R)\ 2 <C(l + \k\ 2 )\\{I-P}u\\l r 

Similarly it holds that 

\A jm (h)\ 2 < C\\v-^{1 - V}h\\\ v \A jm ({l - P}u)| 2 < C\\{1 - P}uf Lr 

Thus, (|3.18j) follows from taking the real part of (|3.19[) and plugging the estimates of I\, I2 
into it. 

Estimate on c u . We claim that for any < 82 < 1, it holds that 

d t ReY, (Bj({I - P}u) I %/) + (1 - S 2 )\k\ 2 \c q \ 2 
j 

< 5 2 \k\ 2 \Vf + % + |fc| 2 )||{I - P}n|| 2 2 + ^^{I - T>}h\\l v (3.22) 

02 t 02 « 

In fact, similarly as before, from the Fourier transform of ()3. 13j) 

d t Bj({I - P}u) + ikjS = Bj(R + h), 

one can get 

d t (Bj({I - P}u) I ik jC ^ + \kj\ 2 \c q \ 2 
= (Bj{R + h) I ik jC ^ + (Bj({I - P}u) I ikjdtc*\ = I3 + h. (3.23) 
1$ is bounded by 

02 

3 



and from the Fourier transform of (|3.9p 

d t S + -ik ■ b q + — V ikjBMl - PVu) = 0, 

j 

I4 is bounded by 

h< 6 ^\k\ 2 m 2 +^\k\ 2 ^2\B J ({i-p}u)\ 2 . 

Notice that similar to Aj m , it holds that 

|Si(^| 2 <C(l + |fc| a )||{I-P}Si||^, 
mh)\ 2 < C\\v^ 2 {l - ?}h\\ 2 \B 3 {{1 - ?}u)\ 2 < C\\{I - P}u\\ 2 Lr 
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Then, (|3.22[) follows from (|3.23j) by taking summation over 1 < i < n, taking the real part 
and then applying the estimates of ^3 and I4. 

Estimate on a u + nc u . We claim that for any < 83 < 1, it holds that 

dtRe^OC I ik m (a q + nc q ^j + (1 - <5 3 )(1 + \k\ 2 )\a q + ndf 

m 

< \k\ 2 \b^\ 2 + Sfc| 2 |c"| 2 + ^\k\ 2 \\{I - P}n|| 2 2 . (3.24) 

In fact, by taking the complex inner product with ikj(a u +nc u ) and then taking summation 
over 1 < j < n, it follows from (|3.2U|) that 

d t J2( b l I %K + rac ")) + + nc "l 2 + Y^~ ik i® I ik i^ + nc ^ 

j j 

= (-2ik jC U I ikjicfi + nc")) + ^(-i%A jm ({I - P}u) \ ikj(aP + nc s )) 

jm 

+ I ikjd t (a? + nc s )). (3.25) 

j 

Using (|3.2ip . one has 

J2(~ikfi I ikji^ + nc")) = ^(kj ^^^ \ a q + nc s ) = |a" + nc £ | 2 . 
j j 

The first two terms on the r.h.s. of (|3.25j) are bounded by 



S 3 \k\V + nc q \ 2 + ^\k\ 2 \S\ 2 + ^\k\ 2 \\{l - P}n|| 2 2 , 
while for the third term, it holds that 

J^(6f I ikjdti^ + ruP)) = J2Q>j I ikj(-ik ■ 6 s )) = |fc • 6 £ | 2 < |fe| 2 |6"| 2 , 
j j 

where we used the Fourier transform of (|3.7I) : 

d t {cfi + nc q )+ik-b q = 0. (3.26) 

Then, putting the above estimates into (|3.25 j) and taking the real part yields (|3.24p . 

Therefore, (|3.17p follows from the proper linear combination of f|3. 18|) . (|3,22p and (|3.24j) 
by taking < <5i,<52)^3 < 1 small enough and also kq > large enough. This completes 
the proof of Lemma 13.11 □ 



Lemma 3.2. It holds that 

n |a c + nc"| 2 

9t y u h + m 2 

for any t > and k S M n . 



\ + A||^/ 2 {I - P}u\\ 2 L 2 < C\\u' l ' 2 {l - P}/i|| 2 2 (3.27) 
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Proof. The Fourier transform of ()3. 1 j) i and (|3. 1 j) *? gives 

k _ ^ 

d t u + i£ ■ ku + i—n 1 {a u + nc u ) ■ £VM = Lu + h. 
\k\ 

Further taking the complex inner product with u and taking the real part yield 

\d t \\uf L2 + Re f («i(a" + nc a ) • ^VM | u) d£ 
* f JR n V \k\ J 

= Re / (Lu | u)di + Re / (h\ u)d£. (3.28) 
For the second term on the l.h.s. of (|3.28j) . from fj3.26[) . one has 

Re / ( iA( « + nc S).fVM |u] d£ = Re / f -^(a e + nc s ) | -ifc • 6 s ) d£ 

7r« V Fr / Jr" \\ k \ J 

=Re L iw^ + ncS) 1 dt ^ + nc "0 de 

- 1 9t|a c + 7ic e | 2 . 



2\k\ 2 

For two terms on the r.h.s. of (|3.28p . one has 



Re f (Lu | u)d$ < — A|| ^ 1 / 2 {I - P} 



and 



Re f (h | = Re / ({I - P}fc | {I - P}u)d£ 

JR™ JR™ 



< ^||^/ 2 {I - P}n||i 2 + C||z,-V2 { i _ 

where P/t = was used. Plugging the above estimates into (|3.28p gives (I3.27p . This 
completes the proof of Lemma 13.21 □ 

Proof of Theorem I3.lt Let k 2 > be a small constant to be determined later. Define 

E(u(t,k)) = ||u(t,A:)||*a + j^la^ + nc q \ 2 + K 2 ReE free (u(t,k)) 

for t > and k £ W 1 , where E free (u(t, k)) is given by (pTTE|) . Notice from (13^61) that 

|£/ ree (u(t, jfe))| < C(|6"| 2 + |c e | 2 + |a s + nc s | 2 ) 



+ ^(|^({I-P}n)| 2 + |^({I-P}2)| 2 ) 



<G(||Pu||^ + ||{I-P}fi||i|) 
<C\\u(t,k)\\l 2 . 
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Therefore, one can choose K2 > small enough such that 

E(u(t, k)) ~ \\u(t, k)f T 2 + + nc q \ 2 (3.29) 

€ \k\ z 

holds. By further letting K2 > be small enough, the linear combination of (|3.27p and 
(|3.17p implies 

I 2 

uu\2 1 1 «|2 



+A|a s + nc"| 2 < C||^- 1 / 2 {I - P}h\\ 2 L 2. 



Notice that 



V2{I _ p} 2 ||2 + ( \ b u\2 + | C S|2\ + | fl fi + nc 



2 



' L ? l + |fc| 2 



U 1 „„M|2 



1 2 



1 



Then, it follows that 



d t E(u(t, k)) + -M^E{u{t, k)) < C\\ V - l l 2 {l - V}h\\l v (3.30) 



which by the Gronwall inequality, implies 



E{u(t,k)) < E{u{0,k))e T W 

+C f e~^ {t " S \\v' 1,2 {l-V}h{s,k)\\ 2 L2 ds (3.31) 
Jo 5 

for any t > and k G l n . 

Now, to prove (|3.3p and (|3.4p . let /i = so that u(t) = e* B uo is the solution to 
the Cauchy problem (|3.ip and hence satisfies the estimate (|3.3ip with h = 0. Write 
/c Q = jfe? 1 ^ 2 • • • k% n . By noticing 

\\d«e m u \\ 2 + \\d°V x A- 1 P e tB uo\\ 2 

= [ \k 2a \ ■ \\u(t,k)\\ 2 L2 dk+ [ \k 2a \ ■ -^\cfi + nc^\ 2 dk 

<C I \k 2a \ \E{u{t,k))\dk, (3.32) 
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then, from (|3.31j) with h = and (|3.29|) . one has 
||^e tB n || 2 + H^V^A-iPoe^noll 2 



X\kj 2 



<C I \k 2a \e *^ \\u5(k)f La dk + C 



I j,2q I a I k I ^ 

Ll e "W*||p^r (fc)||2 2( ifc. (3.33) 



As in [16] . one can further estimate the first term on the r.h.s. of (|3.33p by 



\k 2a \e T W||^ 



idk 



< 



|fc|<l 



/ 2 

7i | ti — 2\ol — ol'\ 



2 dk + 



\k\>i 



e-^\k 2a \ ■ \\^{k)\\ 2 L2 dk 



< C(l + t)—^^^\\d« uo||| 9 + Ce-2*||a> || 2 , (3.34) 

where the Holder and Hausdorff- Young inequalities were used, and similarly for the second 
term on the r.h.s. of (|3.33p . it holds that 



l.2«| A|fc| 2 . 

'e W *||P n (A ; )|| 2 2 dA : < 



n |fc| 



< C(l + t)" 



J 5 



|£ 2q | -^1 



n . n-2(|q-c/|-l) , \ 

« + 2 l|^«o||| +Ce—* t \\d%u f. 



(3.35) 



Thus, (|3.3p follows from f|3.33|) together with (|3.34p and (|3.35p . Moreover, to prove ([37 
notice Po{I — Po}^o = and hence (|3.4p similarly follows only from (|3.33p and (13.341) 
since the second term on the r.h.s. of (|3.33p vanishes. 
Finally, to prove (|3.5p . let uq = so that 



u(t) = f e {t ~ s)B {I-P}h(s)ds 
Jo 



is the solution to the Cauchy problem (|3.ip and hence satisfies the estimate (|3.31|) with 
uq = 0. Then, similar to (|3.32p . one has 



r t 2 t 

dx e {t - s)B {I-P}h(s)ds + d%V x A- l P Q ^"^{I - P}h(s)ds 
Jo Jo 

u - 1 / 2 {I-P}h(s,k)f L 2dkds. 



< C 



k za e i+W 



Therefore, (|3.5p follows in the same way as in (|3.34p . This completes the proof of Theorem 

We conclude this section by extending Theorem 13. II to the exponential time-decay rate 
in the case of T n as mentioned in Remark 13.21 
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Theorem 3.2. Suppose 



PuqoIx = 



(3.36) 



for any £ £ M. n . Let e <B uo be the solution to the Cauchy problem 
d t u + i- V x u - • £ VM = Lu, 



(3.37) 



(3.38) 



'MudS, / = 0, J>0,ieT n ,^ M n , 

u\t=o = Uq, X £ T". 
Then, there are constants C > 0, A > swc/i i^oi 

1 1 e tB 1 1 z, 2 (T" x ) < C e A *II'"o||l 2 (T"xIR™)) 
/or any i > 0. 

Proof. In fact, on the one hand, (|3.17p and (|3.27p with h = in Lemma 13. II and Lemma 
13.21 still hold for any t > and k £ Z n for the solution u(t) = e* B uo to the Cauchy 
problem (|3.37p in the torus case. On the other hand, from (|3.7p - (|3.9p . one has the total 
conservation laws 



T" 



(a u ,b u ,c u )dx 



f (a u °,b u °,, 



c uo )dx, 



which together with the assumption (|3.36p imply 



Thus, it follows that 



which yields 



(a u ,b u ,c u )dx 



i>0,fc=0 



' , , 9 (|&T + IcT + K + ncT) > l&T + |cT + K + nc 
1 + \k\ z 

for any t > and k £ Z n . Therefore, similar to get (|3.30p . it holds that 

d t E(u(t, k)) + XE(u(t, k)) < 0, 
where E(u(t,k)) is still given by (|3.29p . Then, 

E{u{t,k)) < e- xt E{u{0,k)). 
holds for any t > and k £ Z n . Notice that in the torus case, one has 



„«|2 



(3.39) 



1 



E(u(t,k))dk~ I \\u(t,k)\\ z L2 dk+ I —^\a u + nc u \ z dk 



u{t,k)\\j 2 dk+ / |a u +ncTdfc 



\\u(t,k)\\ L2 dk — ||u(t)|| i2 ( T n xR n). 
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Therefore, the integration of (|3.39[) over fc£Z n gives (|3.38p . This completes the proof of 



4 The nonlinear system 
4.1 Energy estimates 

From now on, we devote ourselves to the proof of Proposition 12,11 For that, assume that 
all conditions of Proposition 12.11 hold, particularly 



is supposed to be sufficiently small throughout this section. Let u be the solution to 
the Cauchy problem on the nonlinear VPB system (|2.2p . (|2.3p and (|2.4j) obtained by 
Proposition 15.21 Here, notice from Remark 12.21 that the solution u indeed exists under the 
assumptions of Proposition [2~T1 Therefore, from (|5.5p . 



holds for any t > 0, where £ (u(t)) and T>{u(t)) are defined by (|5.3p and (|5,4p . respectively. 
Thus, one can suppose that the energy functional £(u(t)) is small enough uniformly in 
time. We also remark that those uniform a priori estimates given in Lemma 15.11 will be 
used later on in the proof. 

In this subsection, for some preparations, we are concerned with energy estimates on 
the microscopic part {I — P}u to obtain some Lyapunov-type inequalities. Recall from [8] 
that {I — P}u satisfies 



Theorem 13.21 



□ 



u o\\H N nLlnz 1 




(4.1) 



d t {I - P}u + £ • V,{I - P}u + V x .$ • V ? {I - P}u 
L{I - P}u + r(u, u) + U- V X ${I - P}u 

-{I — P}(£ ■ V x Pu + V x $ • V^Pn - ±£ • V x <S>Pu) 

-P(£ • V a {I - P}u + V«$ • V 5 {I - P}u - ^ • Va,*{I - P}n) 



(4.2) 



and also the following lemma was proved in |29|. IIP] . 



Lemma 4.1. It holds that 



v 



r(u,u)|| L 2 < Cdjz/ 1 7 n|| L 2 ||«|| L 2 + llull^Hz/ 1 7 v|| i 2) 



/or < 7 < 1, and 



H-fTiill^iv < C||n||^jv, 
T(u,v)\\ h n < C(||«|| H Ar||t;|| H jv + |M| H w||u|| H iv). 
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Firstly, corresponding to (|5.6p . it is straightforward to prove 
Lemma 4.2. It holds that 

^||{I - P}u\\ 2 + AHu^I - P}u|| 2 < C\\V x Pu\\ 2 + C£(u(t))V(u(t)), (4.3) 

for any t > 0. 

Proof. The direct zero-order energy estimate on (|4.2|) gives 
~\\{l-P>\\ 2 + X\W 1/2 {l-P}nf 

< I {T(u,u) 1 {I-P}u)dx+ [ (i^ • V^, ({I - P}u) 2 )dx 

- / ({I - P}(£ • V a Pu + V** • V 5 Pn - h • V,.$P U ), {I - P}u)dx, 

which further implies 
1 d 



^ll{I-P}«l| 2 + A||^ { l_ PH 2 

< c^(«(t))p(«W)lk 1/2 {i - p}«|| + c^Rt))lk 1/2 {i - P}n|| 2 

+C||V :c Pu||||i/ 1/2 {I-PH, 

where Lemma [4.1l and Sobolev inequality were used. Then, (|4.3p follows from the Cauchy's 
inequality, smallness of £(u(t)) and the equivalence w(£) ~ ^(0- This completes the proof 
of Lemma l4.2i □ 



Next, we consider the weighted energy estimates on {I — P}u, whose aim is to obtain 
the weighted high-order Lyapunov inequalities later. 

Lemma 4.3. It holds that 

j t \\w 1/2 {l - PH 2 + A|KI - P}u|| 2 < C[l + £{u{t))]V{u{t)). (4.4) 



Furthermore 



l<H<iV 1<M<N 



<CV(u(t)) + C£(u(t)) Y, l«df{I-P}u|| 2 , (4.5) 

|a|+|/3|<JV 



and 



| £ ||^af{I-P}n|| 2 + A £ l«9f{I-P}n|| 2 

I/3|>1 I/3|>1 
|a| + |/3|<JV |a|+|/3|<JV 



<CV{u(t)) + C£{u(t)) l«df{I-P}u|| 2 , (4.6) 

\a\+\/3\<N 

for any t > 0. 
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Proof. The proof of (|4.4p is similar to that of (|4.3p and hence is omitted. To prove (|4.5p . 
let 1 < |a| < TV and then it follows from d2T2]) that 

dt(d°!u) + £ ■ V x (d£u) + ■ V 5 (d» + vd*u 
= Kd%u + Sg'Va.* • ^v 7 ^! + a°T(u, u) - • V{]«, (4.7) 

where the last term on the r.h.s. denotes the commutator. Then, (|4.5p follows from 
multiplying the above equation by wd^u, taking integration over M 3 x M 3 and then using 
integration by parts, Cauchy's inequality, Lemma I4~T1 and Sobolev inequality. 

Finally, we prove (|4.6p . Fix a, (3 with \a\ + \/3\ < N and \{3\ > 1. For simplicity, write 
z = d%dg{I — P}u. Then, after applying d% <9f to (14. 2p . 2 satisfies 

+ £ • V^z + V x $ • V 5 z + 1/2 = I, (4.8) 

where I is denoted by 

I=h+h+h 

with 

Ji = d%Kd£{l - P}u + r(u, u) + • V*${I - P}u), 

7 2 = -a?flf {I - P}(£ • V*Pu + V x * • V 5 Pu - ^ • V^Pu) 

-8£d%P{£ • V X {I - P}u + V a .$ • V 5 {I - P}u - ^ • V a .d>{I - P}u), 

i 3 = -[af , c • v x mi - p}u - [af , i/(e)]^{i - p} u 

-[^V^-V^I-PK 

Here, notice that J2 contains the macroscopic projection P which can absorb both the 
weight and derivative of velocity variable, and ^3 contains all the commutators in which 
the total order of differentiation of {I — P}u is no more than N. Then, similarly for (14. 7p . 
the direct energy estimate of fj4.8|) yields (|4,6p . This also completes the proof of Lemma 

□ 

4.2 Time decay of energy 

In this subsection, we shall prove (|2.26p and (|2,28p in Proposition 12. 11 For this, define two 
temporal functions by 

StHt) = sup (l + s)? +m £(u(s)), m = 0,l, (4.9) 

0<s<t 

where £{u(t)) is defined by (|5.3p . For simplicity, also define two constants depending only 
on initial data by 

e o = IKII#iv + Ikollli, eo,v = ll^oll^iv + Ihollli + lk 1/2- "o|| 2 - (4.10) 
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Proof of (I2T26D and (pT28l) in Proposition [27D We begin with the proof of (01261) . 
Firstly recall that from ()5.5p . one has the energy inequality 

j t S(u(t))+\V(u(t))<0, (4.11) 
where T>(u(t)) is defined by (j5.4h . By comparing (15, 3p with f)5.4|) . it holds that 

£>(u) + ||PH| 2 + HV^A^^otill 2 > X£(u). 
Then, it follows from (14. lip that 

j t £(u(t)) + X£(u(t)) < C||Pm(t)|| 2 + CIlV^A^^o^WII 2 . (4.12) 

Also recall that the solution u to the Cauchy problem (j2.2p - (12,4j) of the nonlinear VPB 
system can be written as the mild form 

u(t) = e m u + f e Bit ~ s '>G(s)ds, (4.13) 
Jo 

where the source term G given by (|2.19p is rewritten as 

G = Gi + G 2 

with 

Gi = r(«,«), G 2 = -V x ^-V^u+^-V x ^u. (4.14) 

By checking PGi = and P0G2 = 0, one can decompose G as 

G = {I-P}G 1 + {I-P}G 2 + P 1 G 2 . (4.15) 

By applying Theorem 13.11 with the spatial dimension n = 3 to 13H and using the 
decomposition (|4.15p . one has 

HPi^f + llv.A^Po^)!! 2 

< C(l + i)~*(||^o|li 1 + ll^oll 2 ) 



3=1 J ° 

/ (1+t- S )-i(||PiG 2 ( S )|| Zl + \\P 1 G 2 (s)\\)ds 
J 



.1 

+C 



(4.16) 



where the first term on the r.h.s. follows from ()3.3p . the second term from (|3.5p and the 
third term from (|3.4p . One can further compute those terms of G\, G 2 in f)4. 16|) by 

||z,-V2{i _ PjdH^ + Wv^il - PjdH 2 



= ||^ 1 / 2 r(u, U )||| 1 + ||^ 1 /2r(n,n)|| 2 

^ r<\\ 1/2 i|2|i 11 2 1 r<\\ 1/2 i|2 n ii2 

< O || 7 u\\ \\u\\ +C\\u ' u\\ sup||u|| L 2 



< C\\v l l 2 u4 2 £{u) 
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from Lemma 14.14 an d 

||^ 1 / 2 {I-P}G 2 ||| 1 + ||^ 1 / 2 {I-P}G 2 || 2 

<^HV^|li2 nLS o(||V ? n|| 2 + ||^ 2 u|| 2 ) 
<C(£(u ) + \W 1/2 u \\ 2 )£(u), 

and 

llPiGyk + HPi^ll < C||V^||^ nL «||«|| < C£(u). 
Here and hereafter, one can use the uniform bounds of £(u(t)) and \\^ 2 u(t)\\ due to 

sup £(u(t)) < £(uq) < Ceo, 

sup \\^ 2 u(t)\\ 2 < C\\^ 2 u \\ 2 + £(u ) < Ce ,„, 

t>0 

where the time integration of (14.111) and (|4.4p was used. Plugging the above inequalities 
into P~T6|) . it follows that 

\\Plu(t)\\ 2 + \\V x A- 1 P oU (t)\\ 2 

i r* a 

< C(l + t)-5e + Ce^ v / (1 + i - s)-3£(t*(a))d« 



+C 



1 - 



(l + t-a)~*£(u(s))ds . 
JO L - 1 

By the definition of £^(t) in (14. 9j) . it further follows that 

\\Piu(t)f + HV.A-^onWII 2 < C(l + «)~3 {e + eo^W + [4?(*)] 2 } , (4-17) 
where we used 

/"* 3 1 1 

/ (l + t-8)-3(l + 8)-*ds<C(l + t)-*, 

Jo 

f f 3 1 1 

/ (l+t-s)-4(l+ S )-2ds<C(l+t)~4. 

JO 

Due to the Gronwall inequality, (|4.12p together with (|4.17p yields 

£(u(t))<£(u Q )e- xt + C f e- x ^- s \\\VMs)\\ 2 + \\V x K l Vou{s)\\ 2 }ds 

Jo 

< C(l + t)~\ {e + e ^\t) + [£^(t)] 2 } , 

which implies 

e®{t) <c{e + e , u £®(t) + [£^(t)?} , 
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for any t > 0. Therefore, as long as both eo and are small enough, one has 



sup£W(t)<Ce , 



(4.18) 



which proves (|2.26p by the definition of Sso (t) in (|4.9p . 

Next, one can modify the above proof of (|2,26p to obtain (|2.28[) under the assumption 
(|2.27p . In fact, under the additional condition (|2.27|) . (|4.16p can be refined as 



where instead of using (|3.3p to estimate the first term on the r.h.s. of (|4.16p . we used (|3,4I) 
since uq does not contain the hyperbolic component, i.e., uq = {I — P}uq. Corresponding 
to obtain (flTTf) . from the definition of s£\t) in (|4"U]) . (|4TT9l) implies 



Pm(t)|| 2 + WV^Pou^f < C(l + t)-§ {eo + e 0|J ,££>(t) + [^(t)] 2 } • (4-20) 



under the smallness condition of eo and eo v - This proves (I23HD by the definition of (t) 
in 

4.3 Time decay of high-order energy 

At this time, in order to complete the proof of Proposition 12. 1\ it suffices to prove (|2.29|) 
in Proposition 12.11 which gives the time-decay rates for the high-order energy. First of 
all, let us obtain some Lyapunov-type inequalities of the high-order energy on the basis of 
Lemma |4.2| Lemma 14.31 and Lemma l5.1l in the following two lemmas. 

Lemma 4.4. There is a high-order energy functional £ h (u(t)) defined by 



IIPi^f + llv.A^Po^)!! 2 

<C(l + trl(\\u \\ 2 Zi + \\u \\ 2 ) 




(4.19) 



Then, from the completely same proof as for (|4.18p . one can obtain 

sup£«(i)<Ce , 



t>o 



£ h (u(t)) ~ ||{I - P l }u(t)\\ 2 HN + ||VJPiu(*)||? 



(4.21) 



such that 



j t £ h (u(t)) + XD{u{t)) < CllV^Pn^ll 2 
holds for any t > 0, where V(u(t)) is defined in ([5. 41) . 



(4.22) 
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Proof. Define 

= ii{i-pm*)ii 2 + E (\\d>(t)\\ 2 + \\d«v x m\\ 2 ) 

l<\a\<N 

+K 3 £ x>i (u(t)) + K 4 £free(u(t)), (4.23) 

where £ x ^(u(t)) and £f ree (u(t)) are given by (|5.9|) and (|5.1U|) . respectively, and constants 
«3 and K4 to be determined later satisfies 

< AS 3 < K4 < 1. (4.24) 
Notice from the definition (|5.10p of £f ree (u(t)) that 

\£free(u(t))\<C £ ( || d«{I - P}u|| 2 + ||^P u|| 2 + || d« V x 6" || 2 ) 
|ce|<AT— 1 

<c||{i-PMt)|| 2 + c Yl Q\%<i)f + \\fi£v**(t)\\ 3 ), 

l<\a\<N 

which together with (|4.23p and (|5.9p imply 

£ h {u(t)) ~ ||{I - PMi)|| 2 + £ (||9>(t)|| 2 + ||fl£V.<IKt)|| a ) + ks^ («(t)) 

l<|a|<JV 



[I - P}«(t)||^ + ||P n(t)|| 2 , (Hf _ 1) + ||V a! P 1 u(t)||^ ( ^_ 1) 

2 1 ||V7 "D „./-AI|2 



~ ||{I - Pi}u(t)|| HJV + HVajPiU^. 



by taking K4 > small enough and also letting K3 > 0. Thus, (|4.2ip holds true. Moreover, 
under the condition (|4.24|) . the linear combination of (|4.3|) . (15. 7p . (|5.8j) and (|5.1ip in the 
special case when 5$ = from the assumption (jl.4|) yield 



^W*)) + \V{u(t)) < C\\V x Pu(t)\\ 2 + C(£ («(*)) + v^Ri)))Z>(u(t)), 

where T>(u(t)) is given by (|5.4p . Therefore, (|4.22p follows from the above inequality and 
smallness of £{u(t)) as mentioned at the beginning of Subsection 14.11 This completes the 
proof of Lemma 14.41 □ 

Lemma 4.5. There is a weighted high-order energy functional £^(u(i)) and a correspond- 
ing dissipation rate T> w (u(t)), which are defined by 

£t{u{t)) ~ ||{I - PiHt)||^ + \\VJ> lU {t)\\ 2 q{H ^ ly (4.25) 

V w (u{t)) ~ ||{I - Pi}u(t)||^iy + ||VxPiu(*)||? a ,„*-i v (4.26) 

such that 

jfiHt)) + XD w {u{t)) < C\\V x Pu(t)f (4.27) 

holds for any t > 0. 
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Proof. Define 

£t(u(t)) = \\w 1 / 2 {I-P}u(t)\\ 2 + £ W wl/2d >^ 

l<\a\<N 

+ £ n^/^^i-Piuf + ^Ki)), 

I/3|>1 
|a|+|/3|<JV 

where £ h (u(t)) is given by (|4.2ip and K5 > is a constant to be determined later. By 
using (I4.2ip . it is straightforward to check 

£*(u(t)) ~ ||{I - PMt)||^ + HV.P^t)!!^ N + K 5 £ h (u(t)) 

,,2 £ V x } 



~ ||{I - PMt)||^ + ||V x Pu(t)||= 3 ^ + ||P u(t)||^ 

~ IKi-PiHt)!!^ +||v x Piu(t)||2 W _! V 

which implies (|4.25|> . The rest is to verify (j4.27H . Actually, by taking K5 > large enough, 
the linear combination of (fO|) . (f4"3j) . P~6j) and P~22j) yields 

^U*(u(t)) + XV w (u(t)) < C\\V x Pu(t)f + C£(u(t))\\{I - P}u(t)\\l N , (4.28) 



2 



where T> w (u(t)) takes the form 

V w {u{t)) = \\w{l - P}n(t)|| 2 + J2 W wd >( 

l<\a\<N 

+ £ l«^{I-P}u(t)|| 2 +P(u(t)). (4.29) 

I/3|>1 
|a|+|/3|<iV 

From the definition (|5.4[) of T>(u(t)), (14.290 implies that (I4.26P holds true. Since for the 
second term on the r.h.s. of (|4.28p it holds that 



£(u(t))\\{I ~ P}<t)f H N < C£(u{t))V w {u{t)), 

i)2 



then (|4.27p follows from (|4.28|) and smallness of £(u(t)) again due to (|4.ip and smallness 
of £{uq). This completes the proof of Lemma 14.51 □ 

Now, we are in a position to prove (|2.29p . For this, as before, define a temporal function 

by 

<«,(*) = sup (1 + sf^£t{u{s)\ (4.30) 

0<s<t 

where < e < 3/4 is an arbitrary constant, and £^(u) is given by (|4.25p . For simplicity, 
also denote a constant t\ depending only on initial data by 

£ i = IKH wjv + I Kill, • (4-31) 
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Notice £^(uq) < Ce\ which will be used later. 

Proof of (|2.29p in Proposition I2.lt We begin with the Lyapunov-type inequality (|4.27[) 
for the weighted high-order energy functional £^(u(t)). Since 

< CV w {u{t)) 

holds by definitions (P5]) - (f06|) of £*{u{t)) and V w (u(t)), (g^7D implies 

^£*(«(t)) + A£*(u(t))<C||V iB Pu(t)|| 2 . 
From the Gronwall inequality, it follows that 

£t(u{t)) < £^(u )e- xt + C f* e- x ^\\V x Pu(s)\\ 2 ds. (4.32) 

J o 

Notice that (|4.32p together with (|2.26j) imply that £^(u(t)) has at least the same time- 
decay rate with the total energy £(u(t)), that is, 

sup(l + *)3£*(u(t)) < Cei, (4.33) 
t>o 

if conditions of ([2.26|) hold and ||iio||#iv is bounded. In what follows, we shall improve 
(|4.33p up to the almost optimal rate in the sense that for any < e < 3/4, there is rj > 
depending only on e such that whenever eg < rj 2 with eo defined in (|4.10p . one has 

sup(l + t) 2 ^£^{u{t)) < Cei, (4.34) 

which implies (|2.29p . In fact, one can obtain a formal time-decay estimate on the first-order 
energy ||V z Pu(t)|| 2 in terms of £%, oo (*) i n the same way as in (14.17P or (|4.20p . Firstly, 
similar to get (|4.16p . by applying Theorem 13. II to (|4.13p with the help of the decomposition 
(|4.15p for the source term G, it follows that 

llV^t)!! 2 <C(1 +i)-3(||u ||| 1 + l|V^o|| 2 ) 

2 r* 

+ VC / (1 + t - srHWu^VAl - P}G? i (a)||| 1 + II^^V^I - P}G j (s)\\ 2 )ds 
~[ Jo 

/ (l + t-s)-m\V ^^(s^z, + \\V x P 1 G 2 (s)\\)ds , 
Jo 

where G\ and G2 are defined in (I4.14p . Next, it is straightforward to check 

\\u- l,2 V x {l - P}G j {t)\\% 1 + \\v- x l*V x {l - P}G j (t)\\ 2 
<C£t{u{t))[£t{u{t))+£{u{t))} 

for j = 1, 2, and 

WVxP^mzr + llV^PiGa^ll < C[£*(u(t))£(u(*))]5. 



+C 



(4.35) 
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Plugging the above inequalities into (|4.35p gives 



\\v x p u (t)r < c(i + t)-2(\\ U0 r Zl + nv^oii' 



+C I (l + t-s)-?£*(u(s))[£*(u(8))+£(u(8))]d8 

-i 2 







+c 



(l + t- s)-*[£*( U (s))£(u(s))}~2d S 



(4.36) 



Notice that from the time integrations of (|4.1ip and (|4.22p as well as the definition (|4,31l) 
of ei, it holds that 

sup[£^u(t))+£(u(t))]<Ce 1 , 



t>o 



and also recall from (l4TT5j) . and (|4TTUD that 

£(u(t))<Ce (l+tyK 
where cq is defined in (|4.10p . Then, it follows from (|4.36p that 



liv^t)!! 2 <c(i + ty 



k + HV^oll 2 ) 



+Cei / (l + t-s)-l£^(u(s))ds 
Jo 



-Ce 



(1 + 1 - s)-i(l + s)-^[£^(u(s))]2ds 



which further from the definition (|4.3U|) of £^ ^ (t) implies 



liv^p^)!! 2 <c(l + ^r 



k+llv^o|| 2 ) 



+Ce 1 £t, 00 {t) f (1 + t - + s)- 2 ^ds 

Jo 

[ (l+t-s)-|(l + s)~3- ( |- e) ^ 

Jo 



n 2 



(4.37) 



+Ce £^ OD (t) 
Notice that for < e < 3/4, 

\l + t - S )-I(l + s )- 2 (!- e )ds < C(l + t)" 2( 3" e ), 
/ (1 + t - s)"t(l + s)-^- e ^ds < C 6 (l + 

■/ 

where C > in the first inequality can be taken uniformly in < e < 3/4, while C 6 > 
in the second inequality has to tend to infinity as e goes to zero. Then, it follows from 
(Q71) that 

||V x Pu(t)|| 2 < (l + tJ-d-^tCeo + CCei + C^^Ct)]. 
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Combined with the above estimate, (|4.32p gives 

S*(u(t)) < {I + t)-^[Ce l + (Ce 1 + 0^)8^)1 
for any t > 0, where max{^(«o), eo} < Ce\ was used. Thus, it follows that 

^ )OO (i)<C ei + (C ei + C e e )^ iOO (t), 

that is, 

since ei is small enough. Therefore, given any < e < 3/4, one can choose n = l/ v / 2C e 
so that whenever €q < i] 2 , it holds that 

£^, 00 (t)<Ce 1 , 

for any t > 0, which proves (|434"j) and hence ([Z32D by the definition (|O0j) of £^ j00 (t). 
This also completes the proof of Proposition 12.11 

5 Appendix 

In this appendix, we shall state the existence of the stationary solution and its nonlinear 
stability for the VPB system (11.1 1> — () 1 .2f) . For that, let us define the weighted norm || • || \y m ,^> 
by 

\\g\\ wr ~ = sup + £ \ d %9(x)\, 

for suitable g = g{x) and for an integer m > and 9 > 0. [8] proved 

Proposition 5.1 (existence of stationary solutions). Let the integer m > and 9 > 0. 
Suppose that \\p — l\\ w m,oa is small enough. Then the following elliptic equation with 
exponential nonlinearity: 

admits a unique solution <p = <j)(x) satisfying 

W^Ww™' 00 — c\\p~- i\\ w ™><*>, 

for some constant C. 

From Proposition 15.11 above, it is straightforward to check that the VPB system (II. 1D - 
(|1.2|) has a stationary solution (/*, <!>*) given by /* = e^M, = 4>. To state the stability 
of the stationary state (/*,<&*), set the perturbation u = u(t,x,^) by 

/ = e*M + VMu. 
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Then u satisfies the perturbed system: 



Att\x\ J R 3 
with given initial data 




(5.1) 



u(0, x, = u (x, = /0 JL , (x, e R 3 x E 3 . 




(5.2) 



Here Lu and T(u,u) are defined by in the same forms as in (|2,5p and (|2.6p . respectively. 
[8] also proved 

Proposition 5.2 (stability of stationary solutions). Let N > 4. Suppose that \\p — 
1||™.jv+i,oo is small enough. Then, there are the equivalent energy functional £(•) and 
the corresponding energy dissipation rate V(-) defined by 



then the Cauchy problem (|5. 1|) - ()5.2j) of the VPB system admits a unique global solution 
u(t, x,£) satisfying f(t, x,£) = e^M + \/M.u(t, x, £) > 0, and 



The following lemma was obtained by [8] in the proof of uniform a priori estimates for 
the global-in-time stability of the stationary solution. 

Lemma 5.1 (a priori estimates). Let all conditions of Proposition WM hold and let u be the 
corresponding solution, and (a u , b u , c u ) be defined in f)2.8|) - ()2 . lQj) . Denote 6<f, = \\(f)\\ w N+i,oa . 
Then, the following uniform a priori estimates hold for any t > 0: 

(i) zero-order: 




—£{u(t))+XD(u{t)) < 0. 



(5.5) 



~(imi 2 +iiv.*ii 2 -2 [ 



e~^\b u \ 2 c u 



dx + A 




M 3 xR 3 



v(£)\{I-P}u\ 2 dxd£ 



< C{5^ + y/E^)))V(u(t)); 



(5.6) 



(it) spatial derivatives: 




< C(5^ + yf£(^tj))V(u(t)y, 



(5.7) 
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(Hi) mixed spatial-velocity derivatives: 

±e Xi t(u(t)) + x J2 II v({)\d«dP{i-p} u \ 2 dxdt 

at J-^ JJr3 xR 3 « 

\a\+\P\<N 

< C(5^ + yf£(aJty))V{u{t)) 

+C \\d%V x (a u ,b u ,c u )f + C Y, |k 1/2 ^{I-PM| 2 , (5.8) 

|a|<JV-l \a\<N 

where 

N 

£ x , i {u{t)) = Y J C N ,k £ ll^{I-PH 2 , (5.9) 

k=l \f3\=k 

\a\ + \(3\<N 

for some proper positive constants CV k> 

(iv) macroscopic dissipation: There is a temporal free energy functional £ f ree (u(t)) in the 
form of 

SfreMt)) =*0 £ £<A;(^{I " P}«), + dtf)) 

\a\<N-l ij 

|ce|<JV-l i 

|a|<JV-l i 

- £ <^(o" + 3c«),^V a! -6 u ) (5.10) 

|a|<JV-l 

for some constant kq > 0, such that one has 

j t £ free (u(t)) + XD mac (u(t)) <C H 5 x{I " P}^f + £(u(t))V(u(t)), (5.11) 

\a\<N 

where T> mac (u(t)) is the macroscopic dissipation rate given by 

V rnac (u{t))= £ \\d«V x (a u + 3c u ,b u , Of + \\a u + 3c u \\ 2 . 

\a\<N-l 
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